Nonadiabatic transitions between adiabatic surfaces: phase diffusion in 
superconducting atomic point contacts 
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Motivated by experiments with current biased superconducting atomic point contacts the general 
problem of nonadiabatic transitions between adiabatic surfaces in presence of strong dissipation is 
studied. For a single channel device the supercurrent is determined by the diffusive motion of the 
superconducting phase difference on two Andreev levels. These surfaces are uncoupled only in the 
adiabatic limit of low to moderate transmissions, while for high transmissions curve crossings are 
important. Starting from a general master equation of the full density matrix an approximate time 
evolution equation for the populations on the adiabatic surfaces in the overdamped limit is derived 
from which the relevant observables can be obtained. Specific results for the case of atomic point 
contacts are in agreement with experimental observations that cannot be explained by conventional 
theory. 

PACS numbers: 74.50.+r,05.40.-a,74.25.Fy, 74.78. Na 



I. INTRODUCTION 

In recent years curve crossing problems regained con- 
siderable attention in contexts such as quantum infor- 
mation processing, superconducting mesoscopic circuits, 
or wave packet dynamics in molecular structures. The 
paradigmatic situation has been formulated by Landau, 
Zener, and Stiickelberg [Ij, i2, iSj] already during the hey- 
day of quantum mechanics: the energy difference be- 
tween the diabatic energy levels of two states, say |1) and 
1 2), changes linearly in time as vt with constant veloc- 
ity V through a crossing region while transitions between 
these states with frequency Aq/H occur. The probabil- 
ity for staying in state |1) in the infinite future when 
starting in |1) in the infinite past is given by the fa- 
mous Landau-Zener-Stiickelberg formula P = exp(— t;) 
where the parameter g — ^/^^/{hv) controls whether the 
dynamics happens to be adiabatically (5 3> 1) or dia- 
batically [g ^ 1). Various modifications of this original 
setting have been studied, in particular, situations where 
the diabatic energy levels vary due to the dynamics of 
an intrinsic degree of freedom. For condensed phase sys- 
tems this degree of freedom typically interacts with a 
dissipative environment so that due to its diffusive dy- 
namics the crossing region (the so-called Landau-Zener 
range) is not traversed ballistically but rather stochasti- 
cally. One prominent example is the charge transfer be- 
tween donor and acceptor states in molecular structures, 
a process which is of relevance for molecular electronics 
Here, in most cases environmental degrees of free- 
dom are fast compared to electronic transitions meaning 
that one is close to the diabatic limit. Important observ- 
ables are then transfer rates between diabatic surfaces 
and time dependent populations on these surfaces. Since 
in general exact analytical solutions are impossible, ap- 
proximate treatments have been developed e.g. based on 
master equations d, H, 0, H] ■ 

Much less attention has been paid to the diffusive dy- 



namics on adiabatic surfaces and the impact of non- 
adiabatic transitions. One realization of this situation 
can be found in mesoscopic physics, namely, in the charge 
transfer through atomic size contacts between super- 
conducting leads 0, [l^, [OL In fact, this device ex- 
tends Josephson weak links [12l | to junctions with variable 
transmissions even close to 1 with the relevant degree of 
freedom being again the phase difference between the su- 
perconducting reservoirs. Theory shows that for an 
isolated contact with one channel the Cooper pair cur- 
rent is carried by two Andreev bound states the energy 
levels of which depend on the phase difference and the 
transmission probability r of the channel. For perfect 
transmission (r = 1) the minimal gap between the levels 
closes, while it tends towards the superconducting gaps 
in the reservoirs for r — + 0. 

For contacts embedded in an electrical circuit the sit- 
uation is more complex. First, the phase interacts with 
the electromagnetic modes of the surrounding and its dy- 
namics becomes diffusive. Second, the contact can be 
externally driven by a current or a voltage bias. The 
charge transfer is then a combination of Cooper pair tun- 
neling, responsible for a supercurrent peak, and multiple 
Andreev reflections (MARs) [H, [H, [3l , responsible for 
quasi-particle transfer leading to a rich subgap structure. 
The first process dominates for low voltages, the second 
one for higher voltages, while in a crossover region a sub- 
tle interplay of both processes appears [T3| ■ 

In this paper we concentrate on the current biased sit- 
uation. Then, the phase diffuses on tilted Andreev levels 
which constitute the adiabatic levels of the system. Due 
to the small capacitance of the contact and the large 
admittance of the circuit the Brownian motion of the 
corresponding fictitious particle is overdamped so that 
this adiabatic approximation is justified. The theoretical 
framework to calculate averages of supercurrent and volt- 
age, respectively, follows that of overdamped Josepshon 
junctions according to Ambegaokar, Halperin and 
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Ivanchenko, Zirbermann [l9|. Ac-driven atomic point 
contacts have been treated in the same way in Ex- 
perimentally, this picture has been verified e.g. by mea- 
suring the supercurrent peak and its maximum (switch- 
ing current) as functions of temperature pll. [2^. Prob- 
lems arise, however, close to the ballistic limit (r close to 
1), when the minimal gap between the adiabatic levels 
becomes comparable to typical diffusion times. Conse- 
quently, nonadiabatic transitions enter the game. Note 
that these transitions can also be interpreted as the pre- 
cursors of MARs. In [2l[ the standard Landau-Zener- 
Stiickelberg formula has been applied to correct for those 
transitions, but failed to capture the experimental data. 
In this paper we provide a consistent theory which ac- 
counts for the overdamped Brownian motion on coupled 
adiabatic surfaces. Such a theory may in turn be of rel- 
evance in other contexts as well and particularly allows 
to calculate steady state currents and diffusive transition 
rates between adiabatic surfaces. For atomic point con- 
tacts its only limitation is that MARs are not included 
explicitly in the theory so that the full crossover between 
supercurrent peak and the subgap structure cannot be 
described. For this case (and a voltage biased contact) 
an effective approach has been outlined recently in p7| . 
which, however, is still somewhat ad-hoc as it combines 
a microscopic description of the MARs without environ- 
mental modes and the phase dynamics in presence of the 
circuit modes. 

The article is organized as follows: In Sec. |TT]we briefly 
describe the standard adiabatic model and introduce the 
relevant notation. A general approach to capture nona- 
diabatic transitions during diffusive motion on adiabatic 
surfaces is formulated in Sees. IIIII and IIVI which is spe- 
cialized to overdamped Brownian motion in |Vl There 
also an explicit expression for transition rates is derived. 
Results for the case of atomic point contacts are discussed 
in IVII The role of quantum fluctuations is addressed in 
Sec.lvni 



II. ADIABATIC PHASE DYNAMICS 

According to the conventional theory of a current bi- 
ased point contact with superconducting leads [ll|, [HI 
transport is described by the phase dynamics on poten- 
tial surfaces corresponding to energy surfaces of Andreev 
bound states. For a single channel with two Andreev 
states the latter ones read 



E±{^) = ±A5^1-rsin2((^/2) 



(1) 



with As being the superconducting gap, r the transmis- 
sion probability through the contact, and (p the phase 
difference between the superconducting reservoirs in the 
leads (see fig. [T|). In the basis of right and left mov- 
ing waves (Bogoliubov-de Gennes theory) the underlying 
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FIG. 1: Energy surfaces E± of the bound Andreev states 
for a single channel with transmission r (solid) together with 
the diabatic surfaces Vi,2 (dotted). The minimal energy gap 
between the Andreev levels at the Landau-Zener point (p = n 
is indicated by the arrow. See text for details. 
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FIG. 2; Equivalent circuit of the resistively and capacitively 
shunted junction (RCSJ) model including a weak link. 



Hamiltonian of this two level system is given by [22 
cos{ip/2) 



Ha 



\/l~^sin(v3/2) \ 
Vl - rsin(v5/2) - cos((^/2) J 



(2) 



so that transitions between "diabatic" surfaces V1/2 = 
± cos((^/2) are mediated by off-diagonal couplings which 
are maximal at (p ~ n corresponding to a minimal energy 
gap 2A5VI — T between the adiabatic surfaces H]). This 
gap closes in the ballistic limit t — 1. The current carried 
by each of the surfaces ^ follows from 



1 dE±{p) ctAs 



sm (p 



tpa dip 



2h 



(3) 

with the reduced fiux quantum ipQ — h/2e. 

A realistic description of transport across the contact 
has also to account for its electromagnetic environment. 
According to the RCSJ model one writes for the circuit 
in fig. [2 



Idc + In = P>oC 



d<2 



^0 



1 dkp 



(4) 
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where an applied bias current consists of a dc-component 
Idc and a fluctuating component /„ which obeys Johnson- 
Nyquist characteristics, i.e., 



(/„(i)) = , (/„(t)/„(0)) 



R 



(5) 



Further, C denotes the capacitance, R the resistance, and 
/(<p) the supercurrent through the contact. This way, the 
phase dynamics is equivalent to the Brownian motion of 
a fictitious particle with mass m = tp^C and friction con- 
stant 7 = 1/RC. In the adiabatic approximation the 
motion of the phase is assumed to be much slower than 
any other relevant process meaning that the supercur- 
rent /((/?) is simply the Boltzmann weighted sum of the 
individual currents /±, i.e.. 



/(^)=/_(^)tanh [PE+{ip)] 



(6) 



with inverse temperature /3 — \/k^T. The effective po- 
tential felt by the fictitious particle is then 



-Idc ifi - 



(7) 



In the actual experimental set-up the capacitance is neg- 
ligible leading to strongly overdamped phase dynamics. 
Hence, the adiabatic approximation is indeed justified as 
long as the adiabatic surfaces are sufficiently separated 
from each other (a more detailed condition will be given 
below) . 

To calculate current and voltage across the contact 
it is much more convenient to work with the probabil- 
ity distribution W{ip^ ip) in phase-space corresponding to 
the Langevin dynamics In the overdamped domain 
the only relevant information is carried by the marginal 
probability distribution n{ip,t) ~ J dipW{ip, tp) and the 
Fokker-Planck equation for the full distribution reduces 
to a Smoluchowski equation of the form 



dn{(p,t) _ R d 
dt Pq dip 



U'{p) + fceT 



d_ 

dip 



nicp,t). (8) 



with U'{p) — dU{ip)/dp. Now, for fixed applied bias 
current I^c the dynamics tends towards a steady state 
distribution Ust {f) for longer times from which the aver- 
age supercurrent through the junction Ij{Idc) = {I{ip))st 
and the mean value of the voltage across the contact 
U{Idc) — 'Po{'p)st can be calculated. Upon varying Idc 
one gains the supercurrent voltage characteristics Ij(U). 



III. DISSIPATIVE TWO LEVEL SYSTEM 

The adiabatic treatment is based on a separation of 
time scales, which may become critical in certain ranges 
of parameter space. The goal of the next two sections is 
to derive formally exact equations of motion for the den- 
sity matrix in the adiabatic basis of the two level system 
thus including nonadiabatic transitions between them. 



This in turn allows to systematically go beyond the sim- 
ple approach presented above. The corresponding anal- 
ysis is completely general and applies to all situations, 
where diffusive dynamics on single adiabatic surfaces 
tends to break down. Accordingly, we slightly change 
notation for the relevant degree of freedom from ip to q 
to stress the close analogy to the dynamics of a fictitious 
particle. 

The starting point is a standard system -I- reservoir 
Hamiltonian of the form H = Hq + Hw + Hb , where 



2m 



V,{q) Ao(q)/2\ 
Ao(g)/2 V2{q) h 



(9) 



describes the dynamics on two coupled diabatic surfaces 
Vi.2 which cross at the Landau-Zener (LZ)-point q* , i.e. 
Vi{q*) = V2{q*), but are separated sufficiently away from 
q* by energies larger than the thermal energy scale k^T. 
Further, 



Hw — 



N 

-q^CnXn + q 

n=l 



N 

2 



(10) 



denotes the interaction of the relevant dynamical degree 
of freedom to a heat bath 



N 



pI 

2m„ 



1 ^ 

n=l 



(11) 



consisting of a large number of harmonic degrees of free- 
dom. In the continuum limit the spectral density of this 
oscillator bath is chosen to be ohmic below a cut-off fre- 
quency LOc, i-e. Ji^^) — mjuj with friction strength 7, and 
zero otherwise. The dynamics of the reduced density 
matrix p{t) = Trsiptotit)} can be represented formally 
exactly in terms of path integrals [l^l . In the domain of 
higher temperatures though, an explicit time evolution 
equation can be derived Jl, [2^ , namely, 

- = [Ho, p] - [q, [q, p]] ~ - [g, {p, p)\ , 

(12) 

where { , } denotes the anti-commutator. The detailed 
conditions for the validity of this approximation are: 
hPujc,fiP7 ^ 1- Note that this regime particularly in- 
cludes the domain of high friction 7/ix'o 3> 1 at suffi- 
ciently elevated temperatures where ujq is a typical sys- 
tem frequency scale. 

For the two level problem the above master equation 
has been the starting point for various approximations 
[2^ 1. In particular, in the high temperature domain 
uJoh/S 1 it is convenient to represent it in terms of 
the Wigner transform [3] 



m.,P,t)./|;.e--'/. 



Pit) 



q 



(13) 
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This gives in leading order in h the Wigner-Fokker-Planck 
equation 

dW _ p dW ^^,dW ^ ( ^ pW 
dt m dq ^ dp ^ \dp^ 

■-mkBT^]+—^{a,,W} 



2h 



dp^ J A dp 



(14) 



where we introduced the sum and the difference of the 
diabatic surfaces, respectively, 

V+ = (Vi + V2)/2 , V-=Vi-V2. (15) 

For practical applications, however, this set of time evolu- 
tion equations for the matrix elements Wij{q,p,t)^i, j = 
1,2 is of limited use only. Namely, it turns out that sin- 
gularities appear for the off-diagonal Fokker-Planck op- 
erators due to highly oscillatory terms V-/h, particularly 
in cases where the so-called reorganisation energy Er is 
large compared to ksT, hcoo Q- Here, Er — \V-^ {qn-iin)\ 
is the energy needed to switch from the minimum g,„in of 
the lower diabatic surface to the higher lying one. In case 
of the two Andreev levels introduced in the previous sec- 
tion, we have Er « Ag. The strategy to proceed is then 
to use approximate solutions for the off-diagonal elements 
in the limit fceT, Hloq Er to derive a set of two coupled 
effective equations for the diagonal densities 0, H, HI] . 
Coupling terms describe the impact of adiabatic transi- 
tions between the diabatic surfaces the strength of which 
is measured by a LZ type of factor AQ/|V"_i|. This way, 
populations and rate constants for the electron transfer 
between donor and acceptor states in molecular struc- 
tures have been calculated numerically and compared to 
analytical findings, see e.g. 0, Hi- The diabatic surfaces 
in this situation are assumed to be harmonic so that the 
system, starting initially from the donor state, relaxes 
towards thermal equilibrium in the long time limit. 

In the present problem the situation is different 
though: here, the relevant physical observables are de- 
termined as expectation values on adiabatic surfaces. To 
calculate them from the diabatic representation is possi- 
ble in principle, but fails in practice since it requires the 
accurate knowledge of the full diabatic density matrix, 
which is not available as discussed above. In particular, 
in a steady state off-diagonal elements typically do not 
vanish. Hence, we follow another route and first trans- 
form the full master equation ([H]) to the adiabatic basis 
before further approximations are applied. 



DYNAMICS OF THE DENSITY MATRIX IN 
THE ADIABATIC BASIS 



IV. 



The master equation (|T4)) describes the dynamics of 
the Wigner distribution in the so-called diabatic basis, 
in which the Hamiltonian becomes diagonal for the arti- 
ficial spin held fixed, that is for vanishing coupling Aq. 



In this section the corresponding master equation in the 
adiabatic basis is derived, the latter one obtained by di- 
agonalizing Hq ^ — {p^ /2m) 1, i.e.. 



Vi Ao/2 
Ao/2 V2 



(16) 



The corresponding unitary transformation is given by 



where (t){q) = (1/2) arctan [Ao/VL(g)]. One finds for 
^'ad ~ UH^U'' the expected result 



H'ad - 



V+ + Vad 

V+- Vad 



with the adiabatic surfaces V-^ ± Vad where 



Vad = \xlV^ 



(18) 



(19) 



describes the gap between them. 

To transform the full Hamiltonian iJo , the unitary op- 
erator must also be applied to the kinetic part. From 
U pU^ = p — hk{q)ay with 



k{q) 



dcp{q) A[,V^_(g)-Aor (q) 



dq 



^Vad{qf 



(20) 



we obtain in the adiabatic representation for Had = 
UHqW 



Ha 



-1 



H'ad 

2m 



h^k{qf 
2m 



^ay{p,k{q)} 
2m 



(21) 



In the adiabatic approximation discussed in the first sec- 
tion, the last two terms are neglected. Here, we retain 
them to find together with U {Hw + Hb)U'^ = Hw ~\- Hb 
the exact form of the master equation (fT2l) for the density 
matrix in the adiabatic basis pad — UpW as 

dPad i,„ , mjkBT 

—Q^ = -j[Uad,Pad\ ^ [qAl^PadW 

Pad}] + —[q,{k{q)ay,pad}] -(22) 

Obviously, the last term accounts for nonadiabatic tran- 
sitions in the adiabatic dynamics. Its role becomes 
much more transparent in the Wigner transformed mas- 
ter equation, where we have 



dW+ 



P 



dt 

dW^^ 



dt 



C+W++ + k{q)^ {W+^ + 
m 

C-W— - k{q)— (W+_ + W-+) 



= CoW+^-k{q)^{W++-W^^) 
dt m 



(23) 
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and likewise for the complex conjugate W |- 

witli the operators 



m oq 



d , 92 
'mkBT-—r 

(24) 

where rj = 0,+,—. Transitions between the adiabatic 
surfaces V++Vad and V+ — Vad occur via transitions to ofF- 
diagonal elements. The corresponding coupling strength 
is given by k{q)p/m which contains both the inverse of a 
LZ type of factor k[q) measuring the differences between 
the diabatic forces compared to the adiabatic energy gap 
and a dynamical factor p/m. Away from the LZ point q* 
one has |V"_| ^ Ag so that the first factor is small and 
together with a slow motion in p nonadiabatic transitions 
are negligible. In the LZ range around V- = 0, however, 
the crucial quantity is |l/i/Ao|, which may become large. 
For instance, for an atomic point contact the LZ range is 
located around (/? = tt where | K!, / Aq | oc 1/ \/l — t. As a 
consequence, the adiabatic approximation breaks down 
for r — > 1 and the dynamics of the density matrix fol- 
lows from the full set of equations only. As in case of 
the diabatic equations (|14p . a direct numerical evalua- 
tion of is in most cases prohibitive though. Namely, 
the off-diagonal elements W± tend to oscillate strongly 
for increasing Vadi i-e. away from the LZ-domain. Phys- 
ically, this reflects the fact that these off-diagonal ele- 
ments are relevant only in a domain around the LZ point 
V- = 0, while outside they are effectively washed out. 
The idea to proceed in the overdamped limit is thus sim- 
ilar as in the diabatic representation: one formally solves 

for the localized dynamics of and W y. and inserts 

these results into the equations for the diagonal elements. 
Eventually one arrives at two coupled effective equations 
of motion for the latter ones which are amenable to nu- 
merical approaches. Corresponding approximations are 
then adapted to the diffusive dynamics on the adiabatic 
surfaces. 



V. OVERDAMPED DYNAMICS 

In the regime of strong friction the dynamics in posi- 
tion is slow, while equilibration in momentum occurs on 
the fast time scale I/7. This separation of time scales 
allows for an explicit elimination of the off-diagonal ele- 
ments in (|23p to obtain a set of equations of motion for 
the populations alone coupled by an effective position 
dependent transition factor. 



A. Population dynamics 

We start by writing ~ exp(2«T4di/ft) W+_ so that 
for the off-diagonal elements in ((23)) one has 



p 2it 



VLW+-. 



(25) 



Now, the propagator of the bare system obeying dQ j dt = 
£oS is to be calculated in the overdamped limit 7/0^0 ^ 1 
with ti>o being a typical system frequency. Accordingly, 
we look for times within the window I/7 < i < 7/(^9, 
where the upper bound follows from the fact that off- 
diagonal elements are determined by the dynamics on a 
much shorter time scale ?i/Ao <C 7/a;Q (see below). This 
way, one finds (see appendix) 



\/2TTmkBT 



1 ^-[p'PQ{t)f/2mkBT ^ 

(26) 

Further, the order of magnitude of the last term in 
(1^5)1 can be estimated in the LZ range as pV^j^t/m ^ 
{{pQ/fn)V^^)p{h/Ao), which can thus be neglected if 



ksT loq 



^0 1 



< 1 



(27) 



Note that for hio^P <^ 1 the above condition is stronger 
than /i/Aq ^ l/^o and thus the relevant one. Now, 
solving (P5|) formally and plugging the result into the 
equations of motion (|23p for the diagonal elements one 
arrives at 



dt 



= L+W++ + 2Ke J dtj dqdpg{q,p;q,p;t-t) 



X ^k{q)k{q) e-2«[^-(?)*-^-(?)fl/?^ 
x[Ty__-W^++](q-,p,i). (28) 



and likewise for W In the overdamped limit the 

Wigner distributions factorize according to 



W±{q,p,t) 



1 



\/2TrmkBT 



exp 



P 



2mkBT 



n±{q,t) 



so that the marginal distributions in position 



(29) 



(30) 



readily follow together with (|26l) from p8)) . We thus 
gain the central result of this work, namely, a time evo- 
lution equation for the adiabatic populations including 
non-adiabatic transitions, i.e.. 



n± {q, t)= J dpW±± {p, q, t) 



9n+ 



= C- 



pp 



k{qf 



pp 



2Vao 



(n_-n+).(31) 



with the Smoluchowski operators 



= CoW+,-^k{q)e^p(^Vadt]iW++-W..) 
ot m \ n 



^ 7717 dq 



(32) 
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and 1] = +,—. Here, Cpp{uj) denotes the Fourier trans- 
form of the momentum-momentum correlation function 



{p{t)p{0)) of a local harmonic oscillator given 



by 



Cpp{t) 



m h 



coth 



cos(ajt) 
(33) 



in the regime of high temperatures where the real part 
dominates [27] . System information is carried by the dy- 
namical susceptibility X = x' + *x"j the imaginary part 
of which becomes in the strong friction limit independent 
of the local oscillator frequency 



/// N 1 7W 



and 



2m? 



(34) 



(35) 



Thus, we see from (|31l) that the nonadiabatic coupling be- 
tween the adiabatic surfaces contains in addition to the 
Landau-Zener factor k{q) dynamical information in terms 
of the power spectrum of the momentum correlations at 
the local transition frequency = 2Vad{q)/fi between the 
adiabatic surfaces. Since the relevant domain for transi- 
tions is the LZ range where the gap between the adiabatic 
surfaces becomes small, one has fio = 2Vad{q*)/fi 7 
and in most cases also iloh(3 < 1 which justifies the high 
temperature analysis given above. A simple extension of 
the result PT|) to lower temperatures {Uhf3 > 1) uses the 
full expression (|33p (see also next section). 



B. Transition rates 

In this section we derive a transparent expression for 
the transition rate between the adiabatic surfaces as de- 
fined by the loss of population in the steady state to the 
other surface if initially only one, say the lower, surface is 
populated. This analysis is independent of the previous 
one and exploits time-dependent perturbation theory. It 
thus applies only if the nonadiabatic coupling between 
the adiabatic surfaces remains sufficiently small. It turns 
out though that this approach provides direct insight into 
the non- adiabatic processes captured in PT|) . 

One writes Had = Hd + Wq with 



H 



H+ 





and the "perturbation" 



p2 fi^k{q)'' 



2m 2m 



iW 
-iW 



(36) 



~^<yy{pM<l)}- (37) 
Zm 



The propagator in the interaction picture is up to second 
order 



Gi{t) 



1 



dt' Wi{t') 



dt" Wi{t')Wi{t").{38) 
The initial density matrix is chosen as 





p_ 



(39) 



with |— ) being the eigenstate of the lower surface and 
P-{q) the corresponding initial distribution in q so that 
the population P-{t) = Tr{|— )(— | pi[t)} is given by 



1- PjRe 



dt' 



' dt"Tr|(GLVK'G+)(i') 
{Glw'G^)[t")-p^] . (40) 



Accordingly, the change in time is obtained to read 

P-{t)~-j^ drTr i^e-'^-^/^We'^+^/^Wp-j , 

(41) 

where one approximates 

exp{iH+T/h) w exp(ii7_r/ft) exp{i2Vadiq)T / h) . (42) 

For sufficiently long times, but short compared to the 
relaxation time, we have P-{t) ~ — F with the formal 
expression for the transition rate 

(43) 

Now, using 

exp{-iH^T/h)W' exp{iH^T/h) w —p(-T)k(q) (44) 

m 



one finds 
F « 



drTr 



kiq) 



p(-T)e*2^°<^^/'"Wp_ 



« 1" dr Tr |p(r)p(0)e-*2K.r/n fcg! | (45^ 

In the overdamped limit the total density matrix p fac- 
torizes in an equilibrium distribution for the momentum 
and a time dependent distribution in position. Hence, 
one arrives with (13311 at 



F= / dr /d<z'Cpp(T)e-^2K-^/'-'^n_(g') 

Vad 



dqk{q) 



2 



coth 



kuT 



- 1 



XX ( ) n_{q). 



(46) 
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As expected, the transition rate contains a transition fac- 
tor appearing also in (|3ip . particularly, the Fourier trans- 
form of the momentum-momentum correlation for strong 
friction. Note however that the full set of equations re- 
mains valid also for larger couplings where the perturba- 
tive approach fails provided ([27|) applies. 



VI. 



RESULTS FOR AN ATOMIC POINT 
CONTACT 



With the general result (pij) at hand we now return to 
the system already addressed in the first section, namely, 
adiabatic dynamics of the phase in single channel super- 
conducting atomic point contacts. We then have for the 
populations on the adiabatic surfaces 



dt 

dn-{(p, t) 
dt 



where the individual diffusion operators read C± — 
{R/iPof3)d^J± with the flux operators 



The phase dependent coupling is 



(48) 



R 



1 



^Vifi [I 



\^/2)Y[l + {2RCVad/hf 



(49) 

and contains Vad — As-^/l — Tsin((p/2)2. 

Before we proceed let us first estimate the range of 
validity of the above equations of motion according to 
the condition (|27p with the minimal Ao((^ = tt) = 
2A5-\/l — For typical experimental parameters [2lLl22| 
'j-'o/t ~ 0.04, and k-gT/ As ^ 0.5 one must obey 1 — t > 
10^'* meaning that the approach only fails for transmis- 
sions extremely close to the full ballistic regime t — 1. 
Note that in this latter limit JC{ip = tt) oc 1/(1 — t) 
diverges. The above restriction can also be understood 
from a different point of view. Namely, the approach dis- 
cussed here has to assume that the instantaneous volt- 
age ifQip produced when the phase diffuses with finite 
velocity through the LZ range is sufficiently smaller than 
the superconducting gap As so that the massive produc- 
tion of quasi-particles is suppressed. The force experi- 
enced by the phase in the LZ domain is V^^ which for 
T — > 1 becomes large in the vicinity of ip — n according 
to {tt — (fi)/ VI — T. Consequently, the width of the LZ 
range tends to vanish for t —^ 1. Thus, estimating the 
typical variance in phase to be on the order of the thermal 
length scale, one derives precisely the above condition. 
Hence, the present formalism is not able to capture the 
full crossover to larger voltages where multiple Andreev 
refiections (MARs) play the dominant role. However, it 
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0.99 



0.998 



T 



FIG. 3: Transition rate according to (146 |l and scaled with 
the frequency fio = Ao(7r)/ft vs. the transmission for various 
inverse temperatures As/ksT = 1 (solid), 2 (dashed), 5 (dot- 
ted). Other parameters are Ec/As ~ 2, 'yh/As ~ 11, and 
/dc<^o/As — 0.3. See text for details. 



does reveal the impact of nonadiabatic transitions in the 
I—U range where the maximum of the supercurrent peak 
is observed. This is discussed in detail below. 

For this purpose, one calculates the mean supercurrent 
through the contact from the steady-state populations 
n^±{4') obeying n^{0) = nj_*(27r) according to 

(Ij) {Idc) = - r VUv) K(^) - nl*(^)] . (50) 
Jo 

The corresponding voltage across the contact is propor- 
tional to the mean phase velocity ([/) ~ <po(<p) and one 
has 



(U) (Idc 



R 



(51) 

The equations (I47|) are now solved numerically with the 
thermal initial conditions n±{t = 0) = exp[— /3(V+ ± 
Vad)]/Z and for fixed values of the external bias cur- 
rent. Throughout the rest of the paper we use typical 
experimental data as in [U . 

We start with the transition rate between the adi- 
abatic surfaces and compare it with the relevant coupling 
frequency at the LZ point JIq — 2AsVT^-T/h. As long 
as r/r2o ^ 1 the adiabatic approximation applies, while 
for r/ilo — 1 nonadiabatic transitions play a dominant 
role. Note that T/Hq diverges for r ^ 1. From fig. [3]one 
observes that in the low temperature range the adiabatic 
approximation is well justified even for transmission very 
close to 1 . The opposite is true for somewhat higher tem- 
peratures, where the ratio T/Hq becomes of order 1 or 
larger. In this domain nonadiabatic transitions must be 
taken into account and the standard Smoluchowski-type 
of phase diffusion on single Andreev levels breaks down. 
This can also be seen from the steady state populations 
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FIG. 4: Steady state populations for the n_ [upper solid 
(blue) line] and the n+ [lower solid (red) line] Andreev lev- 
els for a transmission r = 0.99 according to pip . Also 
shown are the populations according to the dynamics on 
uncoupled surfaces (dashed lines) [2ll . [2^ . Parameters are 
fceT/As = 0.8,<^o/dc/As = 0.3, and ft7/As = 11. 
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FIG. 5: Average supercurrent vs. average voltage in units of 
7?7o with /o = (eA5/fi)(l 
rent of the junction. 



\/l — r) being the critical cur- 
Parameters are the same as in fig. |4l 
Squares (red) denote the results for the coupled dynamics, 
circles (blue) for the standard dynamics on averaged surfaces 
(see text). 



depicted in fig. S) For high transmissions the popula- 
tions around the LZ point differ substantially compared 
to those for isolated dynamics. The supercurrent peak as 
a function of voltage is shown in fig. [5l For comparison 
results of the effective approach outlined in Sec.|TT]are in- 
cluded as well, where a mean potential surface consisting 
of a thermal average of the two Andreev levels is used in 
a standard Smoluchowski equation [see ([7|) , ([8])] . In accor- 
dance with the above discussion non-adiabatic transitions 
do increase the supercurrent Ij compared to the case of 
the dynamics on averaged surfaces, because more (less) 
population is residing on the upper (lower) level in the 
down-hill direction. Experimentally, the variance of the 




keT/As 



FIG. 6: Maximal supercurrent as a function of temperature. 
Parameters are chosen as in fig. |4] Squares (red) depict the 
data for coupled dynamics, circles (blue) for the standard 
dynamics on averaged surfaces (see text). 



switching current (essentially the maximal supercurrent 
-^max) with temperature has been found to deviate sub- 
stantially for high transmittin g ch annels from the stan- 
dard Smoluchowski prediction [2ll [2^ . In fact, with in- 
creasing temperature /max has been observed to be larger 
than the predicted values, for very low temperatures 
smaller. The simple strategy to include non-adiabatic 
transitions by using the standard Landau-Zener formula 
mentioned in the Introduction, however, failed to capture 
this effect consistently. As illustrated in fig.[S]the coupled 
dynamics does indeed describe the observed phenomena 
qualitatively. The increase in /max for rising tempera- 
tures can be attributed to the larger transition rates seen 
in fig. El A quantitative comparison with the experimen- 
tal data necessitates a more careful analysis of the actual 
circuit (see also below) and will be presented elsewhere. 



VII. ROLE OF QUANTUM FLUCTUATIONS 

The theory developed so far is basically a classical one. 
As discussed at the end of Sec. IVl quantum effects in the 
momentum-momentum correlation can be effectively ac- 
counted for by working with the full correlation (j33p resp. 
its Fourier transform. This seems not to be consistent in 
a strict sense since quantum fluctuations in position are 
completely ignored. Here, we will provide arguments why 
this strategy is justified and to what extent the latter ones 
may appear in the theory. 

The classical Smoluchowski theory requires not only 
strong friction 7/0^0 ^ 1 but also 7/1/? ^ 1. The op- 
posite limit where the quantum scale for friction by far 
exceeds the thermal energy scale, i.e. 7/1/? 1, has 
been studied in [1, [H, [29[. The essence is this: In 
the so-called quantum Smoluchowski regime to leading 
order the classical Smoluchowski equation still applies. 
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Quantum fluctuations appear in form of a modified dif- 
fusion coefficient ksT -> ksT/ll - AV"{q)/kBT] wfiere 
A = (?i/m7r7) ln(7ft/3/27r) for ^hp > 1. Typically, A is 
small and describes deviations of equilibrium fluctuations 
in position from its classical value {q^) (3 — {q^) i3,c\- In con- 
trast to this squeezing in position, fluctuations in momen- 
tum are large and fully quantum mechanical in agreement 
with the uncertainty principle. It is thus justified even 
for ^h[3 1 to use at least to leading order the classi- 
cal Smoluchowski equation for the dynamics in position, 
but the quantum version of the momentum-momentum 
correlation. To next order, quantum fluctuations in po- 
sition are accounted for by the above replacement of the 
diffusion coefficient in ((3T|) . 

One may wonder in which regime, classical or quan- 
tum, superconducting point contacts are operated. For 
this purpose one realizes that ^hj3 = fiEc /irp = n where 
Ec = 2e^/C is the charging energy and p = R/Rq with 
= h/ie^ js^l- In the overdamped limit one always 
has p ^ 1. For a contact with a capacitance in the fF 
range and R of the order of 20017 one then gets k « 
40/T[K] so that even for temperatures of a few K, we 
have K S> 1. The conclusion is that the phase dynamics 
of the junction is taking place in the quantum Smolu- 
chowski regime. We note that quantum fluctuations cor- 
respond here to charging effects and display Coulomb 
blockade physics [3Q] . Question is why this has not been 
observed yet. There are two answers. The first one is 
based on the above argument that even in this regime to 
leading order everything is classical. The second one is 
based on an analysis of the real circuitry. There, addi- 
tional capacitances Cg of typically a few pF are placed 
in parallel to the weak link meaning that it sees a more 
complex admittance with additional voltage fluctuations. 
These latter ones are classical since hp/RCg <C 1 and 
may mask the quantum fluctuations. If either or both of 
these arguments apply needs a more careful study of the 
experimental situation and goes beyond the scope of the 
present work. Anyway, atomic point contacts may be 
ideal test-beds to tune an overdamped system between 
its classical and quantum regimes. 



the maximal supercurrent is qualitatively explained. The 
classical theory may be extended to include quantum 
fluctuations. 
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Appendix 

Here, we give a brief account on how to obtain the 
propagator for the off-diagonal elements in in the 
overdamped limit. In this domain and in the time win- 
dow 1/7 ^ t <C 7/wq {ujq is a typical frequency of the 
system) the dynamics of the propagator Q{q,p; q,p; t) ac- 
cording to dQ I dt = CqQ is determined by only local prop- 
erties of the potential V+. Thus, we may use a local 
harmonic approximation V^{q)' ~ y^{q) + yl{q){q — q)- 
With ^ (9)/to2 and / = ^^{g) the operator £0 in 
takes the form 



-^|- + (/ + mf72g)|- 
m aq op 



+7 



, 32 
p + mkBT—r 
op^ 



(52) 



The propagator for this harmonic problem with the ini- 
tial condition Q{q,p\q,p,0) = 6{q — q)5{p — p) is well- 
known 



1 



2Tr\/ FG - m 



exp 



GR^ - 2HRQ + FQ'^ 

2{FG-m) 



(53) 



where the full time dependence is captured by the corre- 
lation functions 



VIII. SUMMARY 

We developed a consistent approach to describe in the 
regime of strong friction non-adiabatic transitions be- 
tween adiabatic surfaces. This results in a set of cou- 
pled equations of motion for the adiabatic populations 
(I3ip . which may be of use in a broad range of contexts. 
An explicit expression for the transition rate reveals that 
the strength of these transition is controlled by the in- 
verse of a LZ-type of factor measuring the diabatic forces 
at the LZ-point relative to the energy gap between the 
adiabatic levels and the momentum-momentum correla- 
tion. Explicit results are obtained for high transmitting 
superconducting atomic point contacts. In particular, 
the experimentally observed temperature dependence of 



F 

G 
H 



and 



{l-e-^* [2(7Viy2)ginh2(,,t/2) 



+ {-f/iy)sinh{iyt) + 1]} , 
mfceT {1 - e-''* [2{-f^/i^^) smh^{iyt/2) 
-(7/i.)sinh(z/t)-Kl]} , 

(47A:Br/i/2)e-T* sinh2(^V2) , (54) 



R = q -qR.it) , Q=p-PQit). 



(55) 



Here, v = — 4122 and the latter two functions carry 
the mean dynamics of position and momentum, i.e, 



= e-^*/2 
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PQit) 



pe 



vt 



X sinh 



sinh 



vt 

2{f + ^1X1^1^] 



(56) 



Now, within the time window I/7 ^ t <C 7/^^^ the cor- 
relation H becomes of order I/7, while the product FG 
is larger than by a factor 7/: ^ 1. Accordingly, one 



shows that the q and the p-dependence in Q factorize. 
Further, for the small parameter e = t/{'-f/fl'^) <^ 1 the q 
dependence becomes a Gaussian sharply peaked around 
q = q with a width of order e thus representing effectively 
a 6{q — q) contribution. In contrast, the p — p fluctua- 
tions are of order 1. This way, one arrives at the result 
specified in (|26p . As expected the dynamics in position 
is frozen within the time window, while the momentum 
equilibrates to its instantaneous value around pq (t) . 
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